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Abstract. The spin glass behavior near zero temperature is a complicated matter. To get an easier access to 
the spin glass order parameter Q(x) and, at the same time, keep track of Qat, its matrix aspect, and hence 
of the Hessian controlling stability, we investigate an expansion of the replicated free energy functional 
around its "spherical" approximation. This expansion is obtained by introducing a constraint-field and a 
(double) Legendre Transform expressed in terms of spin correlators and constraint-field correlators. The 
spherical approximation has the spin fluctuations treated with a global constraint and the expansion of 
the Legendre Transformed functional brings them closer and closer to the Ising local constraint. In this 
paper we examine the first contribution of the systematic corrections to the spherical starting point. 

PACS. 75.10.Nr Spin-glass and other random models - 64.70.Pf Glass transitions 

1 Introduction replica trick, that is all the thermodynamical information 

is encoded in the n ^ limit of the disorder-averaged 



The infinite range Sherrington-Kirkpatrick (SK) model is 
defined by [T]: 

l.Af 



replicated partition function: 



^ ^ ^" = Tr,„ exp 1^-/3 ^HsK(Ta)j (2) 

'HsK = -i^^JijcncTj (1) 

ij where /? = 1 /T is the inverse of temperature and, as usual. 



where the (Ti's are ±1 Ising spins and the couplings Jtj (• • •) the average over the quenched disorder. Perform- 
independent Gaussian random variables with null mean ing this average, and introducing the auxiliary symmetric 
and variance equal to 1/N . This model is solved with the replica matrix Qab = j; X^iLi '^ia o'ib with a ^ b, leads to 
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ID: 



ab 



27r 



(3) 



with the effective Lagrangian (density) [2]: 

ab 

Q[Q] = In Trcr„ exp f — ^ Qab ctq 

^ ab 



(4) 
(5) 



The normahzation factor in eq. ^ gives a sub-leading 
contribution for iV — > oo and is omitted in the following. 

In the thermodynamical limit, iV — + oo, the value 
of the integral in ^ is given by the saddle-point value 
—n(3fsK = ^[Q], where Qab is evaluated from the station- 
arity condition: 

5 



SQal 



■C[Q] =0, a^b 



that is, 



iab = {<7a<7b) 



(6) 



(7) 



The average is taken with respect to the weight 

If instead, one had in mind to write the effective La- 
grangian governing a short-range {e.g. nearest- neighbor) 
system, one would consider: 



I n dQabi^) exp ( - ^ E E(^'' + 1) Qabip) 

i,ah p ah 

+ ^ Q[Q{i)] + \Y. ^-bii) Qabi^)] ■ (9) 



Here Qab{p) is the (space) Fourier Transform of the site- 
dependent replica overlap matrix Qab{i), with i the site- 
index, Hab{i) is an external (unphysical) field that cou- 
ples to Qab{i), and n[Q{i)] is given for each site i by dS]). 



Properties of yy[i?], and of its Lcgendre Transform r[Q], 
defined via: 



w[H] + r[Q] = IJ2 ^-b^^ 



(10) 



and, 



SQab{i) 

are obtained by writing: 



r[Q] = Hab{i) 



(11) 



Qab{i) = Qab + SQab{i) 



(12) 



with Qab given as the (site-independent) saddle-point va- 
lue of £, eq. ([6]). An expansion of ([9]) in powers of SQab{i) 
shows that no linear terms {i.e. tadpoles) remains by vir- 
tue of ([5]). The quadratic term is given by: 

^^5]5g.6(p)M'^''^^''(p)'5gcrf(p) (13) 

p abed 

where M is the Hessian matrix whose eigenvalues are the 
so-called bare masses of the SQ correlation functions {i.e. 
the inverse propagator). With these notations, one has: 

6^ 



M^'^^''{p) = {p' + l)S^l,^-T' 



-f2[Q] . (14) 



SQab^Qcd 

The higher order terms are in turn the couplings (cubic, 
quartic,...) of the fields SQ. 

Whatever the adopted viewpoint, one cannot avoid the 
difficult construction of the i7 functional as given in ^ . 
Parisi and others j5l6l7| have shown how to obtain solu- 
tions with R steps of Replica Symmetry Breaking (RSB) 
and in particular with i? — > cx), and how to construct 
equations satisfied by Q{x), the continuous limit of the 
order parameter Qab for i? — > oo 8 . These equations can 
be solved in the full low temperature phase [9]. Clearly, 
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however, for very low temperatures or for null tempera- endre Transform and the stability of the saddle-point. In 
ture, the problem is delicate, as one can experience when Sections [3] and [4] we discuss respectively the (trivial) ze- 
expanding C(''') in powers of [see eq. ([S])]. The prop- roth order approximation and the first non-trivial approx- 
erties of the Parisi solution Q{x) at very low temperatures imation, the two-replica approximation, of the functional 
have attracted some work in last years. The analysis has f2. To make the paper also accessible to uninitiated read- 
been carried on either directly on the continuous limit, i.e. ers, we have deferred to appendices, besides a couple of 
taking first i? — > oo an then T ^ 1 jlO|llj . or on the T ^ detailed calculations, a reminder concerning some of the 
limit of large i?-step replica symmetry breaking solutions technical tools used in the text (Replica Fourier Transform 
[12113114115] . In both approaches, however, properties of and double Legendre Transform), 
the Hessian remain elusive since once one is working with 

Q{x) the matrix structure of the order parameter Qab is 2 General formalism 

lost. On the other hand taking first T ^ and then i? » 1 

2.1 Equation of motion: 

one encounters problems following the delicate commuta- 

tivity of the two limits, extensively discussed in Ref. [16, . spherical models, the spin variables a, are continuous: 

These two faced difiiculties on the road to obtaining the t^'^^^ "^^^ «Pi^^ becomes an n-dimensional Gaus- 

Hessian (and eigenvalues) have motivated this work. The ^^^^ '^-^tegv^^. Here, the Ising spins are discrete variables 

starting point is the observation that, if the strict "focal" ^ overcome the discrete nature of the spin, one 

constraint cr^ = 1 is replaced, as a first step, by the global introduces a constraint field. 

n 

constraint — n {a kind of spherical approximation TrCT(- • •) = (• ■ •) 

in replica space) then, the temperature behavior gets sim- 
pler while keeping the matrix structure of the problem. 

~°° a=l a=l 

local constraint introducing higher and higher correlations 
generated under a double Legendre Transform on the 12 
functional. In the resulting approximated functional for The auxiliary variables Xa control the fiuctuations of the 
n, the null temperature limit becomes non-singular, and, continuous spin variables around the Ising values ±1 
at the same time, one keeps the matrix structure of Qab, '^^P'^^^ ^^^^ ^^^''''g^ as T -> 0. 

leaving access to the Hessian inserting (fT5| into the expression ^ , the functional 

C[Q] can be written as: 

The outline of the paper is as follows: in Section [2] we 2 i.n 2 

C[Q] = -^y Ql^^ + f2^^[Q]-!^ (16) 
describe the general formalism based on the double Leg- 4 4 



a=l aa=±l 



^ /' + 00 

n / [S{<Ja + l) + S{aa-l)]{---) 



X exp 



(•••)• (15) 
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where 



no+[Q] = lim 



= lim In 



+ CXD 



dX 

{2mY 



-Slrr.X] 



(17) 



and 5[cr, A] is given by, 



S[CT,X] 



ah 



(18) 



The constant term — follows from the (definition 
of the) diagonal terms Qaa = 1 of the overlap matrix. 
We recall that when introduced, the overlap matrix is de- 
fined only for a 7^ 5, the reason being that for Ising spins 
cr^ = 1. The diagonal terms Qaa are then usually (tac- 
itly) set to zero, as in the previous Sections. This choice 
may not be the most suitable for spherical models where 
sums over spin values are replaced by Gaussian integrals. 
We thus (re)define the overlap matrix with the diagonal 
terms equal to 1. The constant term comes out then for 
consistency with eqs. (|4]) and ([5]). 

The functional fi^ [Q] can be expressed in terms of the 
field expectation values and two-point correlation func- 
tions by introducing linear and bilinear auxiliary fields 
coupled with cr and A, and performing a double Legendre 
Transform with respect to these fields |17I18I19] . See also 
Appendix [D1 In absence of an external field coupled to a a 
the average {ua) vanishes. Moreover, one can check that 
(caAb) = [5D]. Therefore the relevant field averages and 
correlation functions are [21] 



and 

{aaCTb) = Gab, (XaXb) = Afe • (20) 

All averages are evaluated with the weight exp(— ^[cr. A]). 

After the double Legendre Transform has been taken, 
and the auxiliary linear and bilinear fields set to zero, the 
functional i7efQl comes out to be 1181: 



f2^[Q] = -S[0,X] - ^TrlnG^i - ^TilnA-^ (21) 
-^TtD-^G -^TtA-^A- IC2[G,A] + n 



where, 

{D-\b = XSab-P^Qab, {A-')ab = -eSab ■ (22) 

Here, /C2[G, yl] is given by the sum of all the two-particle 
irreducible (2-PI) vacuum graphs in a theory with vertices 



1 



a \ 



(23) 



and propagators 



Gab — 



Aab = 'VWWVWW . (24) 

a b 



a b 

The first terms of the diagrammatic expansion of JC2 [G, A] 
are: 



l^2[G,r]= ^y\/\/\/\/\rh 




(25) 



a b 

By construction, when the auxiliary fields are set to zero 

the double Legendre Transform is stationary with respect 
to the variations of A, G and A. This leads respectively to: 



(Aa> - A 



(19) 



Gaa — 1 



(26) 
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{G-^),^ = {D-^),^ + 2-^K2[G,A] 



(27) 



2.2 Stability 



d 



{A-^)ai> = 2—-JC2[G,A] 

OAab 



Equations (P7| and are known as the Dyson equation 
for the two-point functions G and A with self-energies: 



(28) To decide whether the solutions of the equation of mo- 
tion are stable or not, one has to expand the effective 
Lagrangian to second order in 5Q: 



Kb = -2 



Kb = -2 



d 



dGab 

d 
dAai 



K.2[G,A] 
/C2[G,yl] 



(29) 
(30) 



C{Qab + 5Qab{p)}^ C'^^'^Qab} 

+ ^Y.^Q-bmGab-Qab) 



ah 



E E 5Qab{p)M''''^''''{p)5QM (33) 

p abed 



These Dyson equations result from the Legendre Trans- '^i*'^ Hessian matrix: 



form, and give exact relations for G and A. By solving 
them, the propagators Gab and Aab become functionals of 
Qab that, substituted back into (HU, give Q[Q] as func- 
tional of Qab alone. By using now the saddle-point equa- 
tion ^ one can write the self-consistent equation for Qab, 
that reads 



Qab — Gab, a^b . 



(31) 



With the choice Qaa — 1 eqs. (|26)) and (|3T|) can be merged 
into the single equation of motion 



Qab = Gab Va, h . 



(32) 



We stress that, while in the following we use the station- 
arity conditions in this form, one should keep in mind that 
the diagonal and off-diagonal terms of these equations fol- 
low from two distinct conditions. The diagonal terms fol- 
low from the double Legendre Transform, and our choice 
Qaa = 1- This is an exact relation. The off-diagonal terms 
follow from the saddle point eq. (O and hence eq. (|3ip is 
valid only in the thermodynamic limit. 



^-^-'^(p) = (p2 + l)^Kr 

Let us introduce: 

p^ab,cd 
ya 



SGab 
SQcd 



SGab 

SQcd 

b.cd _ SAab 



SQcd 



which through the equations of motion ([?7|) and 



(34) 

(35) 
(36) 
sat- 



isfy 



E '^aeGfjf + G^jG^^ 
ef L 

52 



E 



SGabSAef 



E 



SGabSGef 

s^ 



IC2 



X 



ef,cd 



"ab,cd 
rp2 



(37) 



SAabSAef 



IC2 



y 



-E 



/Co X'^f'"'^ = 



^SAabSGef 

Solving ([55)1 for y, and inserting it into (|57|) . one has 



(38) 



E 

ef 



C i 1 T 1 

^ae ^bf 



E 



— 1 z' ' — 1 

af ^be 



SGabSAgh 

gh " 13 

"ab.cd 
y2 ■ 



SGabSGef 

S^ 

'HAjSGef 



IC2 



■IC2 



p^ef,cd 



(39) 



^ Please note that in eqs. (p7)l - p9)) we used G^^^ for (G ^)ab 
to bypass cumbersome writing. 
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The matrix T) is the one multiplying 3^ in ([38| . At this dinary products of RFT matrices (like in ordinary Fourier 
point, we can now use the equality Gab = Qab, resulting Transform), we thus get [24] : 



from the stationarity condition. 

To solve the self-consistent equations of motion, we 



have to specify the form of the matrix Qab that takes '^J^ 



Ml} = + l)-X[r (42) 
1 

-='!"■ _ _ , 

GkGi 



kl 



(43) 



into account possible breaking of the permutation sym- -^fcf ^ ^ ^ ^kl ^ ^^^-^ 

metry of replica pairs. A standard parameterization has n 

where r = 0, i? and fc, ; = r + 1, i? + 1. By /Cvp-- 

we mean the double RFT on the Replicon component of 

SG \sG with respect the lower indices (cross overlaps) 

u and V. The "hat" denotes the RFT (see Appendix E|. 



been introduced by Parisi 5 6 . The overlap matrix Qab for 
R allowed breaking in the replica permutation symmetry 
consists in dividing Qab into successive boxes of decreas- 
ing size pr, with po = n and pr+i — 1, along the diagonal 
and assigning the elements Qab of the overlap matrix so 
that 



3 Zero order approximation 



Q = Q = Q r — ■ ■ ■ R + 1 (40) ^^'^ zeroth order approximation is obtained by neglecting 

the 2-PI contributions of JC2[G,A]. Combining eqs. (|^^ . 

with 1 = Qr+1 > Qr > ■ ■ Qi > Qo- The notation aCib ~ 

(P7)) and (1221) one is led to 

r means that a and b belong to the same box of size Pr but 

to two distinct boxes of size Pr+i < Pr- The element Qq {Q-^)ab^XSab-l3^Qab- (45) 



must vanish in absence of external fields that break the 
up/down symmetry [22] . In the following, if not explicitly 
stated, this will be always assumed. 

In this formalism the Replicon component of the Hes- 
sian (that usually harbors the most dangerous, i.e., the 



To solve this equation the parameterization of the Qab 
matrix is needed. By using Parisi's parameterization one 
gets for the off-diagonal terms 

{Q-'),^-p''Qr, r = 0,...,R. (46) 



lowest eigenvalues), becomes I23I24I , . . , itti r r, .^ . ,• , 

' ' As shown m Appendix ICl tor any R, this equation only 

M'i^;^ = MZ, (41) admits the solution [55] 

where anb = cnd = ris the common overlap while Qr = Q = I — T, Vr (47) 

u = max{a He, and} and v — max{6 H c, 6 n c?} the cross- 



i.e., the Replica Symmetric (RS) solution. The value of A 
is set by the diagonal term a = 6 of eq. (I45p and reads 



overlaps. Replicon matrices enjoy the property of being 
diagonalized under a double Replica Fourier Transform 

2 

(RTF), and matrix products of such matrices becoming or- ^^^"^ 
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The stability of the solution is ruled by the Replicon eigen- with the matrices Gab and Aab solution of the Dyson equa- 

value [see eqs. (|42])-(|44l)], tions 

7W99(p = 0) = 1-/32 (1-Q)2^0, (49) [G-^b ^ XSab - f3^Qab + GabAab (52) 

(yl-i)„b = ^Gab (53) 

where r = and the (Fourier transformed) cross-overlap 2 

k = I = 1. The solution is then marginally stable, as and Gaa = 1- Stationarity of £[(5] with respect to the vari- 

known for the 2-spin spherical model. ations of Qab leads to the self-consistent equation Qab — 

Let us note a technical point. In this approximation Gab, the form of which is dictated by the ansatz used for 

Aab = —e^ab: SO that the term Trln/l is singular as e ^ the matrix Qab- 

0+. The functional C[Q] is nevertheless regular in this The Replicon component for the two- replica approxi- 

limit. The reason is that neglecting the JC2[G,A\ contri- mation reads: 
bution is equivalent to remove the Ising-like constraint on 



kl 

the spin variables. This is also equivalent to discard the r 1 n -1 



last two terms of S[(t, A] in eq. (jl8p . A proper normaliza- 



Ar - QrAkAiQr 



(54) 



MYfip) =p^ + l 

1 

QkQi 

tion is now needed to regularize the e ^ 0+ hmit of [Q] . with r = 0, ... ,R and k,l = r + 1, ... ,R+ 1. 
This introduces an additional constant term in C[Q] that In the next subsection we examine the existence and 
exactly balances the singularity coming from Tr In A. stability of three possible scenarios, corresponding to three 

particular parameterizations of the Parisi matrix: the Re- 
. ^ ,. . . plica Symmetric, the one step Replica Symmetry Breaking 

4 Two-replica approximation 

and the full Replica Symmetry Breaking ansatz. 

The two-replica approximation consists in taking only the 

first 2-PI diagram in the diagrammatic expansion of /C2[G,yl],|_i jhe Replica Symmetric (RS) solution 

namely. 



In the Replica Symmetric solution all replicas are treated 
^2[G, yl] = #yv\/w\/^ on an equal footing, and the matrices Qab and Aab are 



U-lX '^T.^abAab. (50) 



ab 



The functional C[Q] then reads 

-^[Q] = + jTrlnG+ ^Trlnyl 



given by 

Qab = (1 - Q) Sab + Q, Aab = (^1 " ^o) Sab + Aq . (55) 

Then from eq. (j53p and the relation Qab = Gab we have 
for n ^ 
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that, inserted back into ([52|l . leads to the RS equation 4.2 The one step Replica Symmetry Breaking (IRSB) 

solution 



1-2Q ^ ^2 o 1~2Q' 



2 



A - + 2- (57) 



(1-Q)2 " (1-Q2)2 

^ ^3 In the previous section we have seen that the first cor- 

- ^ — • (58) 

rection coming from IC2 [G, A] is not enough to completely 



(i-g)2 ^ - ■ (i-g2) 



These equations can be easily solved for A and P as func- destabilize the RS solution that exists for T < 1, yet its 
tion of Q e [0, 1], and one ends up with domain of vahdity is reduced to the temperature range 

up to 0.34 .... In this section we investigate the existence 

2 Q 

A = Y~~Q2 ^^^^ of solutions parameterized by the so-called one step RSB 

= ^ + y . (60) ansatz: 



(1-Q2)2 • 



{Q2-Qi)5ab + iQi-Qo)eab + Qo (63) 



By varying Q between and 1 one obtains the solution, 
if it exists, in the whole range of physical values of Q. where the matrix Sab is equal to 1 if a and b are in anyone 
We note in particular that Q vanishes as T ^ 1^, while of the diagonal boxes of size pi = m, and otherwise. 
Q — > 1 as T goes to 0: The diagonal element Qr+i = 1 of the overlap matrix 

writes Q2 = 1. The matrices Gab and Aab have a similar 
^ ^ ^^^^ ^ ^ ■ (^1) structure. 



The stability of the RS solution is ruled by the Replicon 
eigenvalue 



In the absence of an external field that breaks the 
up/down symmetry, Qq = 0. This in turn implies that 
all 0- indexed quantities, such as Aq and Gq, must also 



7\/l99(p = 0) — 1 ^ vanish. It is not difficult to show that Qo = ^0 = Go = 



2 



(i-g)-2 + ^-Q2(^^_^^) 

ACp- is a solution of the self-consistent equations of motion. 

^ (62) 



- 2g - 1 



Inserting the IRSB parameterization into eq. (|53p . us- 
The Replicon eigenvalue is positive at T = 0, where its i^g the condition Gab = Qab, and taking the n ^ limit, 



value is 2, and remains positive up to the temperature have 
T* ~ 0.34..., where the denominator vanishes. For larger 
temperatures the Replicon is negative, and vanishes as 



(64) 



T 1^. Therefore we conclude that the RS solutions 

(l-g^)(l-0^ + mQ2) 



^1 = -T^^^IUT^^T^r—-?^, ■ (65) 



exists for all temperatures below T — 1 , as in the trivial Note that for m ^ {] these relations reduce to those found 
case. But this solution is only stable for T < T*. This is for the RS case. This is not unexpected since if Qo = the 
a rather unusual feature of this approximation. systems breaks down into n/m subsystems of dimension 



A. Crisanti et al.: Low temperature spin glass fluctuations: expanding around a spherical approximation 9 

771, each one having a RS structure, Qi playing the role of It turns out that the IRSB solution appears at the Grit- 
ty . For m —^ we recover the RS solution. ical temperature T — 1 and exists down to the lower crit- 
The other Dyson equation, as of eq. leads to ical temperature Tirsb — 0.656485 . . . where Qi reaches 
,„ , ^ the maximum allowed value Qi = 1. 

(1 + Oi)(l - Qi + m Qi) The stability of the IRSB solution is now ruled by two 

l + 2Q, + {m-l)Ql 

^ {l + Qi)il + Qi-mQi){l + Qj-mQl)'^ ' Replicon eigenvalues. 

As for the RS solution by varying Qi in the range [0, 1] Mfl{p = 0) = 1 - - Qi)^ (70) 
we get, for any value of m e [0, 1] a IRSB solution, if it 

and 

exists. 

To fix the value of m we use the stationarity equa- -^WiP = 0) = 1 (71) 



tion {d/dm)C[Q] = (yielding the so-called static solu- 
tion |26I27I28I29| ). The computation is achieved using the 
formulae of Appendices |^ and [Bl leading to 



1 



1-Qi 



9 



In 



1 - Qi + mQi 
A2-A1 



(l-Ql)-2 + yli-Q2(^2-^)-2 ■ 

At T = 1 both eigenvalues vanish. These two eigenvalues 
cross each other at Tq = 0.6760. Above Tq, Mf^ is the 
lowest eigenvalues and is negative. Below To, A^jj is also 
negative. Therefore the solution is always unstable. 



A2- Ai+mAi 
1 Qi ,1 Al 



ml-Qi+mQi m A2 - Ai + m Ai 43 f^^,, Repij^a Symmetry Breaking (full-RSB) 

-\q\Ai={). (68) . 

^ solution 

By eliminating Ai, A2 and (3 in favor of Qi and m with 

We consider now the solution, if any, with an infinite num- 

the help of eqs. ([64|)-([65ll and ([67ll . we end up with the 

ber of replica symmetry breaking steps (full-RSB). In this 

following equation 

limit the matrices Qab, Gab and Aab are described by func- 
2(1 — Qi+mQi) tions of a single parameter x varying between and 1. To 



1-Ql+mQl 



-mgi(2-Kmgi) = (69) 



find the self-consistent equation of motion of the full-RSB 
solution let us consider the case of R replica symmetry 
breaking steps. The full-RSB solution is obtained as the 



that gives m = m{Qi) for the IRSB case. Solving this limit R 00. 

equation shows that m(Qi) is a monotonically increasing By assuming a Parisi's R replica symmetry breaking 

function of Qi that vanishes for Qi ~ 0, i.e., T — 1, while structure for the matrices Qab, Gab and Aab, the Dyson 

HmQj^i m{Qi) = 41n2 — 2 < 1. equations ([5^ and ([55]) for the non-diagonal terms be- 



10 
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(72) 

r = 0,...,i?. (73) 



These equations are solved through the RepUca Fourier 
Transform (RFT) (see Appendix [X] for proper definition 
and Appendix [Bl for an example). The RFT of eq. (|73p 

reads 

1 1 

(74) 



where we have introduced the shorthand qat — (Qab)'^- 
Inverting the RFT yields in the continuous limit R oo 



Aix) 



ds --^A{s) 



where, 



qix) = 1 



s ds 



Q{x,f+ f ' dss^Qis)^ 
Jx ds 



(75) 



and 



Xc= hm PR, Q{xc) = hm Qr. 

R — >oo _R — >-oo 



(76) 



(77) 



Finally, from eqs. ([7^ . ((75|) and ([75]) we have the equation 
of motion for Q{x) for < x < Xc'- 



ds 3^ = f3' Q{x) + 4 Q{x) r ds ^^4^ (78) 



Q{sy 



where the "dot" denotes the derivative of the function 
with respect to its argument. As done for the IRSB so- 
lution we assumed that all 0- indexed quantities, such as 
Qo = Q{0) etc., vanish. 

The solution of this integral equation is obtained by 
successive differentiation with respect to x. The first yields 



Q{xY 



^^ + 4/ ,,Q(jffi)^4-Q(^)^ 



q{sY 



q{sY 



(79) 



By setting x = we get 

Q{0)^l~Q{x,)+ [ ^ dssQ{s)^T . (80) 
As Q{x) — Q{xc) for X > Xc, equation ([50)) reads 

1 - / dxQ{x) = T . (81) 

"'0 

Notice that the same manipulations applied to the internal 
energy (per site) u, yield 

i ^1 - ^ dxQ^{x)^ = -uT . (82) 

The sum rules (|8T|) and ([82]) are known to be exact for the 
SK model. They are verified here to the first order of the 
2-PI expansion of K,2 ■ 

To solve the equation of motion (|79p we take one more 
X derivative to get 

Q{xfq{,x) 



X = Q Q(x) 



Taking advantage of 



q{xf -^Q{xfQ{xf 



(83) 



Q{Q) = 1-Qc+ / dQ'x{Q') 

JQ 

= l-Qc~viQ) 



(84) 



/•Qc 

q{Q) = l-Ql + 2 dQ'Q'x{Q') 

jQ 

= l-Ql~2Qv{Q) + 2r[Q) (85) 

where Qc — Q{xc) is the plateau value, it is easy to check 
that eqs. ([83)) . ([84]) and (|85p are equivalent to the coupled 
differential equations 

d 



dQ 
d 

dQ 



r{Q)^v{Q) 
v{Q)^x{Q) 



with the boundary condition 



r(Qc) = f (Qc) = 0. 



(86) 



(87) 
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leading to a zero mode: 



11 



Q(x) 




(90) 



Fig. 1. Q{x) versus x for the fuU-RSB solution. The solution 
is obtained with Qc = 0.5 and corresponds to temperature 
T = 0.8076 . . . 

Solving these equations for < Qc < 1, and fixing the 
temperature through eq. (|80p , we obtain the complete so- 
lution Q{x). In Fig. [T] we show the form of Q{x). As found 
for the IRSB solution, the temperature does not vanishes 
as Qc 1, but reaches the finite value limg^^^i T((3c) = 
Tfuii-RSB = 0.7829 . . .. 

At variance with the IRSB solution the continuous so- 
lution is marginally stable since the lowest Replicon eigen- 
values vanishes. Indeed, looking at stability, we have for 
the Replicon eigenvalues: 



1 




1 




iQkQi ^ 





(88) 



QkQl . 



Before letting R oo the lowest eigenvalue is M^J^ ^ . 

r+lr+l 

In the continuum A^52L^ is correspondingly given by, 

X + X + 



(89) 



1 



1 



This is a check of the existence of Goldstone zero modes, 
accompanying the breakdown of a continuous invariance 
group, i.e. the existence of Ward-Takahashi identities |30j . 
From ([55)1 one can also obtain the size of the band on top 



of the zero- modes M'SSL^ by evaluating: 

X+X + 



M 



00^ 

XnXn 



= M 

0+0+ 



XX 



1 - 



1 



1 



( 1 ]^ 

Vl-Qi/ 

1-/32(1-Qi)2 



(91) 



As a result at this level of approximation, one keeps the 
same structure as for the SK model near Tc'. a Replicon 
broad band of modes, bordered below by zero modes. 

5 Conclusion 

In this work, we have introduced, for disordered Ising-like 
systems, an expansion around their spherical approxima- 
tion in replica space via the use of constraint-fields and a 
double Legendre Transform of the free energy functional. 
Successive approximations are given by 2-PI graphs con- 
structed with replicated spin correlators and replicated 
constraints correlators. A striking feature of this approach 
is that one can keep track of the the matrix structure of 
the order parameter (and thus of the Hessian) at all tem- 
peratures. 

We have applied the method to the SK model and ana- 



On ([59)1 . one identifies, from the equation of motion ([75]) . 
that the quantity in the square bracket is equal to T^, lyzed in detail the lowest nontrivial contribution, studying 
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Stable 










IRSB 










Unstable 


PM 


RS 


RS 








1 Stable 


Unstable 








i 











Fig. 2. Phase diagram in the two-replica approximation. The 
RS solution {R = 0) exists for all temperature below T — 1, 
but it is stable below T* = 0.34 . . . only. The IRSB solution 
{R = 1) exists in the temperature range Tirsb = 0.656 . . . < 
T < 1, but it is unstable. Finally the full-RSB solution 
{R oo) exists in the temperature range Ttuii-RSB < T < 1, 
Tfuii-RSB = 0.783 . . ., and it is marginally stable. Above T — 1 
only the paramagnetic (PM) solution with Qab = Sab exists. 
All solutions coincide at T = 1. 

in each case (RS, IRSB, full-RSB) successively, the exis- 
tence of solutions and their stability as summarized in 
Fig. [21 At high enough T, the constraint-field fluctuations 
are sufficient to destabilize the RS solution giving rise, in- 
stead, to a marginally stable continuous solution. At low 
enough T the constraint-field becomes ineffective: a stable 
RS solution is found. 

The ineffectiveness of the constraint-field seems linked 
to the fact that, at its lowest order the 2-PI functional 
is linear in the constraint-field. This leads actually to a 
trivial relationship between the constraint-field and the 
order parameter and a correction to the free energy that is 



linear in temperature, a very peculiar feature of the lowest 
order approximation. This remark motivates future work 
to take into account more generic features of the model as 
in the four replica approximation. 

One of us (CD.) would like to thank E. Brezin and H. Orland 
for useful discussions. A.C. would like to thank the IPhT of 
CEA, where part of this work was done, for hospitality and 
support. 



A Inversion of a Parisi Matrix with the RFT 
formalism 

In the R steps RSB scheme an n x n symmetric matrix 
Qab is parameterized by two sets of numbers. The first set 
(po,Pi, • ■ • ,Pr.,Pr+i), with Pa = n and p^+i = 1, gives 
the size of the boxes the matrix is successively divided 
into. By definition for any finite n the pr are decreas- 
ing integer numbers pr > Pr+i, however in the n — > 
limit they becomes real numbers in the range [0, 1] and 
its order gets reversed, that is, Pr+i > Pr- The second 
set {Qo, Qi, . . . , Qr, Qr+i), gives the value of Qab in each 
box, 



?afc — Qanb—r 



r = Q,---,R+l 



(92) 



The notation a n b = r means that the indexes a and b 
belong to the same box of size pr but to two distinct boxes 
of size Pr+i- The term Qr+i is the value of diagonal term 

Qaa ■ 
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The Replica Fourier Transform Qab of the matrix Qab integral and we have, for Q < x < Xc= pr^ 
is defined by [24131] : 



B+1 ' ' Jo Q(s)2 Q(0)2 

Qfc = ^ Pr(Qr-gr-l), k ^ 0, 1 , . . . , R + 1 (93) 

r=k where 



(99) 



where we have taken Q-i = 0. In the following all quan- 

-.u-A . u A A. Q{x)=Q{l)-Q{x,)+ r dssQ{s) . (100) 

titles with index out ot the allowed range are assumed to J x 

be equal to zero. From the above definition, the useful ^^^^ Q(l) — Qr+i Similarly 
relation 

Qk - 4+1 = {Qk - Qfe-i) (94) ^^"'^^^^ ^ 0(i)-Q(xc) + ^^"'^^"^^ ^^^^^ 

1 f ^ Q(s) Q(0) 



allows the inversion of the RFT Q{1) — Q{xc) Jq Q{s)'^ Qi^)^ 

1 _ 

Or = V — (Ofc-Qfe+i), r = 0,l,...,i?+l . (95) 

— Vh 

'==0 B Computing Tr In Q 

It can be shown that the Replica Fourier Transform 

has the properties of the standard Fourier Transform. In The quantity TrlnQ can be computed by using the iden- 

particular the RFT of the inverse matrix {Q~^)ab is the tity 

inverse of the RFT of Qab' -TrlnQ = lim — ( -TrQ™ - 1 ) (102) 

n m^Q m \n J 

Q-^i^ — ^ . (96) and the results of the previous Appendix. 

Qk 

The quantity Q™ for integer m is the matrix product 

Then, by virtue of eqs. ([M|) and ([55]) we have for r ~ 

of 771 matrices, i.e., a convolution in matrix indexes. As a 

0,---,i?, 



consequence 



EQ k ~ Q A; — 1 / \ 

— ^ ^ — (97) Then, by using eq. (PSj). we have 



(103) 



k=Q Qk+l Qk 

1 1 



Qr+1 



and hence 



t^o ;^TrlnQ=)^- 

fc=0 ^ 

These expression remains valid also in the n ^ limit. 1 f Qk 



In Qfc - In Qk+i 



Pk \ IVlil 

In the limit of i? ^ oo the difference Qr — Qr-i with fe=i 



\n{QR+i-QR) + J2^ In 



r < R goes to zero, the sum can then be replaced by an ' Qi (105) 
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This expression is valid for any finite R. To perform the where 

R ^ oo limit we use the relation ((94|) and expand the ^ rQ fj^Q' 



F{Q) =P^Q- f 
Jo 



logarithm in the sum for smaU Q,. - Q^_i, i.e., J^, dQ"a;(Q") 

ln(S^]=ln(l+p,9l^^] (106) 

\Qk+i J \ Qk+i J ^ 



(112) 



Ik 



r=0 
R 



so that eq. p05)) becomes - ^ [Vr+i ~ Pr)S{Q - Qr) SQr- (113) 



r=0 



\ Q Q 

-TrlnQ = ln((5fl+i — Qji) + "S^ — — By requiring the stationarity of C with respect to varia- 

n Qk 



fc=i 

R tions of Qr and p,. one gets, respectively 



+^ + J20{{Qk^Qk+if) . (107) 
fc=i 

In tlie limit i? — )■ oo the last sum vanishes and we end up 
with 



F(g,)=0, r = 0,...,i?, (114) 
/ dqF{Q)=0, r=l,...,i?. (115) 



^TrlnQ — ln((5(l) Q(x )) T^i^G function F{Q) is continuous, thus Eq. (jllSp implies 

/•a;c h('r\ n(n\ that between any two successive Qr there must be at least 

+ / dx^ + ^ (108) 



n 



q(:e) g(o) 



C Proof of eq. ([47]) 



two extrema of F{Q). If we denote these by Q*, then the 
extremal condition F'{Q*) = implies that 
1 

dqx{q)=T. (116) 



If we neglect the term K,2[G,A\, and use Gab — Qab the 

r. 1 /• 1 l A^T i ^^ i The right hand side of this equation is a constant and 

lunctionai L can be written as [27133] : ^ 

2 hence there is only one solution (or none). Thus the only 

-£ = ^ V Qoh + -Tr In Q + const 

n 2n ^ n solution of eq. ([113 is Qr-i ^ Qr = Q for aU r. This 

= Z?' dQQx(Q) + ^-JR in turn implies that xiQ') = 9{Q' - Q), so that eq. HHj) 

' Jo Jo J^dQ'xiQ') 

+ ln(l-Qii) + const (109) 

where 

R 

X{Q) =P0 + ^(Pr+l - Pr) 0{Q - Qr) ■ (HO) 



reduces to 

/3^Q-(^=0 (117) 
that is, Q = 1 - T. 



The saddle point equations are obtained by varying the D Double Legendre Transform 

above functional with respect to x{Q): 



In this section we give a short summary of the double 



^ ~ lo '^'^^(Q) '^a;((3) (111) Legendre Transform to help to understand the origin of 
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eq. ([4]) . To keep the notation as simple as possible we shall 
consider the case of a single scalar field. The extension to 
more complicated cases is straightforward. 
We are interested in the evaluation of 



In Z = In 



(118) 



for some S[(l)]. To this end we define the generating func- 
tion W[J,K] as [T8IT9] 



,WIJ,K] 



-S[4>]+J4>+^<t,K<p 



Clearly 



InZ = WIJ.K] 



.J=K=0 

From the definition of W\J, K] it follows 



(119) 



(120) 



W[J,K] = {4>)=v 



d_ 
dJ 

±^W[J,K] = ]^{ct>4>)^\{v^ + G) 



(121) 
(122) 



where the average (• • •) is taken with the full weight 



exp 



- s\4>] + -4>K4> 



(123) 



and G is the connected two-point correlator. We can then 
define the double Legendre Transform FYp, G] of W[J, K] 



as 



W[J,K]+r[ip,G] = Jip+'^ipKip+'^KG (124) 



with 



—r[^, G] = j + K^, 



(125) 
(126) 



Comparison of eq. (I120|) with cqs. ()124p - (ll26p shows that 



InZ = -r[ip,G] 



(127) 



where (p and G are solution of the equations 
d 



dip 
d_ 
dG 



r[^,G] = 



(128) 
(129) 



i.e., the value of InZ is equal to (minus) the value of 
r[ip, G] at its stationarity point. 

Up to now we have just used some general proper- 
ties of the Legendre Transform. The usefulness of this ap- 
proach arises from the fact that one has an explicit form 
of r[(p, G]. Indeed one has 



r[ip, G] = S[ip] + iTrlnG"^ + ^TiD-^G 
+/C2b,G]-iTrl 



where 



92 



(130) 



(131) 



and IC2[ip,G] is given by the sum of all 2-PI vacuum di- 
agrams of a theory with propagators G and interaction 
vertices determined by the potential Vi^t given by 



-(f>D-^[ip]<f>+V,^,[(f>;ip] . 

(132) 



This procedure corresponds to a dressed loop expansion 
with vertices that depend on (p and can thus exhibit non- 
perturbative effects even for a small number of dressed 
loops. The stationarity condition on r[ip, G] yields a set 
of coupled (nonlinear) equations for ip and G. 
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